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Abstract 

Many results that are difficult can be found more easily by using a generalization in the complex 
plane of Einstein's addition law of parallel velocities. Such a generalization is a natural way to add 
quantities that are limited to bounded values. We show how this generalization directly provides 
phase factors such as the Wigner angle in special relativity and how this generalization is connected 
in the simplest case with the composition of 2 x 2 S matrices. 
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I. INTRODUCTION 



In special relativity the composition law of parallel velocities appears to be the natural 
addition law for quantities whose values are limited to the closed interval [—1,1], where we 
have set the speed of light c = 1. It is natural to generalize the composition law of parallel 
velocities to the complex plane as 

A = A 2 ©A 1 = A±^L, (1) 
1 + A 2 At 

where A\ and A 2 are complex quantities and where the denominator appears as a normal- 
ization term (if not otherwise stated, A denotes the complex conjugation operation). The 
physical meaning of this composition law is similar to that of the composition law of parallel 
velocities in special relativity. Equation (0Q) shows that no matter what real values we give 
to A\ = V\ and A 2 = v 2 , subject only to v\ < c and v 2 < c, the value of the resulting velocity 
A = w cannot exceed the speed of light c = 1. In the same way, no matter the values of the 
complex quantities A\ and A 2 (subject only to \A\\ < 1 and |A 2 | < 1), the modulus of the 
resulting quantity A cannot exceed unity. 

Because it avoids infinities, such a generalization of Einstein's composition law of veloc- 
ities appears to be a natural addition law in a closed interval. As expected, it reduces to 
the usual addition of arithmetics when the quantities are small. As shown in Refs. [1IJ2IJ3I, 
the use of this composition law quickly leads to important theoretical results and provides 
useful algorithms for computer calculations. The use of Eq. p]) also leads to results which 
converge more rapidly than by using transfer matrices.- 

II. SOME SIMPLE EXAMPLES 

We consider three examplesi^ where the use of Eq. ([1]) is useful. The examples are the 
composition of two non-parallel velocities in special relativity, the reflection coefficient of a 
Fabry-Perot in optics, and the characteristics of a polarizer resulting from the association 
of two successive non-perfect polarizers. 

The composition law for two parallel velocities v± and v 2 in special relativity is (c = 1) 

v i + v 2 /nS 

w = v 2 ®v 1 = -± 2 -. 2 

1 + V 1 v 2 
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The calculation of the resulting velocity of two parallel velocities is straightforward. However, 
it is not when the two velocities are not parallel for which calculations may be tedious. They 
become simple when we consider Eq. ([T]) which is the generalization in the complex plane 
of Eq. ([2]). As explained in Ref. Q, we replace each velocity tj| by the complex number 

^ = tanh|e"% (3) 

where the rapidity czj is related to the modulus of vl by tanh = v j and where the phase 
a, gives the orientation of V{ with respect to an arbitrary axis of the reference frame of the 
observer in the plane of V\ and v~2. The modulus and phase a of the velocity Tu resulting 
from the relativistic composition of V\ and v% is directly obtained 1 by using Eq. ([I]) 

tanh — e iai + tanh — e iQ2 

W = tanh-e ia = V 2 ®V 1 = ^- \ . (4) 

2 i + tanh^e— 2 tanh ^ e JQl 

2 2 

The modulus and the phase of Eq. (J4]) gives respectively the magnitude of the resulting 
velocity Tu (because w = tanh a) and specifies the direction a of ~w in the plane (vi, V2). 

In optics the overall reflection coefficient of a Fabry-Perot interferometer can be obtained 
by taking into account all virtual paths of light inside the interferometer.- The total prob- 
ability amplitude for light to be reflected by the system can also be directly obtained (for 
any number of interfaces) by using Eq. (TjQ). Here 

Ri = r< (5) 

is the complex reflection coefficient of an incident wave on the interface i, where is the 
Fresnel coefficient of that interface and fa is the phase shift corresponding to the propagation 
of light through the same homogeneous layer between two successive interfaces. For two 
interfaces the reflection coefficient of the whole system can be obtained directly by using the 
law dU 

r , _|_ r „ P ic t>2 

R = re** = R 2 © R, = 1 1 \ 2 .. . (6) 

Again, the modulus and the phase of Eq. ([6]) give the overall reflection coefficient and phase 
of the reflected wave. 

Similarly, we can consider the composition of two non-perfect polarizers Pi and P2.— The 
polarizer P resulting from the combination of polarizers Pi and P 2 (in that order) can also 



3 



be obtained from Eq. ([I]). As explained in Ref. y, each polarizer is characterized by 

Pi = tanh — e iai 
2 

where Tj gives the quality of the polarizer (typically, ji = Y^z where Tj is the differential 
absorption rate of the polarizer and z the distance traveled by the light wave inside the 
polarizer; the case of a perfect polarizer corresponds to 7, — > +00), and where a, gives the 
orientation of the polarizer axis with respect to an arbitrary reference axis. The polarizer's 
orientations and the reference axis are coplanar. The characteristics of the resulting polarizer 
P are given by^ 

tanh — e iai + tanh — e ia2 

P = tanhle ia = P 2 ®Pi = ; ^7^- (?) 

2 1 + tanh — e~ tai tanh — e ma 

2 2 

By using the composition law © we easily extract the 7 factor which is the quality of the 
resulting polarizer and its direction a. 

The use of the composition law (pQ) is general and can be applied to any number of coplanar 
velocities in special relativity, to any number of interfaces for the case of multilayers, and to 
any number of successive polarizers. In such cases we have to iterate Eq. (OQ) as relation^ 1 ^ 

A = A n ©(A n _ 1 ©---(A 2 ©A 1 )). (8) 

The successive iteration of Eq. (0Q) yields the desired result. Equation §8$) leads to algorithms 
which are useful for many problems. It is easy to compute A 2 © A\ and then to compose 
the result with A3 and so on. 

As explained in Refs. Handfl the expression for A in Eq. ([8]) can be written down directly 
by using a complex generalization of the elementary symmetric functions of the variables 
Ai, A 2 ,- ■ ■ , A n which are extensively used in the theory of polynomials.- 1 ^ 

Our aim in this paper is to show how Eq. (OQ) is related to 2 x 2 matrices and how it 
provides a simple way to calculate the four elements of scattering matrices (^-matrices). We 
also show how the use of Eq. (OQ) leads naturally to a particular phase, which for the case of 
the special relativity is related to the Thomas precession. 

III. MATRIX REPRESENTATION 

We now explain how the composition law ([TJ) is related to 2 x 2 matrices. 
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FIG. 1: Schematic representation of the linear relations between input and output quantities. 



A. Mathematical definitions 



Consider a physical system (see Fig. [Q in which two physical quantities (the inputs) E^ 
and E^ are linearly related to two other physical quantities (the outputs), E^ ut and E~ ut . 
These relations can be written as a 2 x 2 matrix as 



(9) 



This general representation is for example used to describe a birefringent system with the 





help of Jones matrices (see e.g. Ref. |3|), or to estimate properties of a multilayer stack with 
the Abeles matrices.— It is possible to define using the 2x2 matrix in Eq. ([9]), hereafter 
called [M], the four coefficients TZ + , 1Z~, T + , and T~ (this notation is chosen in analogy to 
the classical reflection and transmission coefficients of a multilayer device) 

Et _n A'' 



T- 



T- 



E, 

El 



out 
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We introduce the variable defined by 



If we use Eqs. ( TTU|) and ffTT]) . it is easy to verify that 

r+ t- = e, 



(10a) 
(10b) 
(10c) 
(lOd) 

(11) 
(12) 



which constitutes a generalization of the Stokes relation which is well known in the optics of 



multilayer devices (see e.g. Ref. Il2|). We also introduce the conjugation operation (denoted 
by the bar): 

72 + = -TT = ^, (13a) 

and 

W = -72+. (13b) 

The conjugation operation does not necessarily correspond to the usual complex conjugation 
(compare Eqs. (j 10 aft and (I13ap ). For hermitian matrices the correspondence does hold. 
With these definitions the matrix [M] can be written as 

T+ \n+ e J 

This form of [M] will be useful in the following derivations. 
B. Composition laws of the 1Z and T variables 

We now focus on the properties of the four coefficients 72+!, 72^, ^21 > an d %i of a 
system characterized by its [M21] matrix, resulting from the composition of two subsystems 
characterized by the two [Mi] and [M 2 ] matrices defined by 

[MJ = -L ( 1 *H and [M 2 ] = -L f ' + (15) 

The [M 2 i] matrix is the result of the product of two matrices: [M 21 ] = [M 2 ] [Mi]. Equa- 
tions f[TTl) - f[T3"j) allow us to express the composition laws for 72^, T^, 72^, and r T 2 \ as 

72+ = 72+ © 72+ = gj^±gt , ( i 6a) 



1 + 72+72 
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^" = 72 2 - © 72r = ^ r t_^ 2 Ql , (16b) 
1 + 72 x 72^ 

T 2 + = r 2 + ©T 1 + = Tl+ ^ 2 1 + , (16c) 
1 + 72+ 72 2 

T 21 = Tf © Tf = Tx 3 . (16d) 
1 + 72 x 72 2 

Equation (fT3|) can be used to show that the denominators in Eq. (TT6T) are the same, 1 + 
72+^ = 1 + 72772 2 . 
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C. Composition law of the 6 variables 



Although the variable has been introduced in an ad hoc way in Eq. fllip . it is interesting 
to find its composition law. We consider two processes characterized by the two variables 
0i and 2 . If we start from the generalized Stokes relation (TT21) and use Eq. ( TTB|) . we find 

e 2 i = T+ T 2l - TZ+ K 2l (17a) 

r + r+ rf t 2 - - [njQ 2 + nj] [n^ + n 2 e 1 ] 

,-.2 ( 17b ) 

1 + ^+^2 

/From Eq. we know that T x + Tf = ©i + TZfTZ^ and T 2 + T 2 = <3) 2 + 71% Tl 2 . Conse- 
quently Eq. (117bj) becomes 

1 + K+K 2 

This expression can be considered as the composition law for 0! and ©2- In Sec. |IV]we will 
give the meaning of for various physical contexts. 

D. S'-matrix 

By definition, the four coefficients 7Z + , T + , 1Z~, and T _ are the four elements of the 
S'-matrix associated with scattering, 




5=1 , . (19) 



Equation ffl6|) shows that the composition of two S matrices can be written as 

S = S 2 o S 1 = 2 1 2 1 . (20) 
V T+ ® T+ K 2 ®K{) 

The use of the composition laws © and <g> give the elements of the S'-matrix without resorting 
to the usual transfer matrices. 



IV. THE G PHASE FACTOR 

We now consider conservative systems described by the unitary matrix \U\. In this context 
the variables are modulus one complex numbers of the form e 1 ^. Our aim is to show that 



the phases associated with the physical modes E + and E~ in Eq. (Q can be written as a 
sum of phases when the two modes are not coupled, plus a phase which is simply expressed 
with the help of the © law. 



A. The composition law in the case of unitary matrices 

The general expression of a 2 x 2 unitary matrix is 

(cos A e lu — sin Xe lv \ 
. e %v , (21) 
sin A e w cos A e m I 

where <p, A, u, and v are real numbers. The overall phase <p can be omitted without loss 
of generality and hereafter we set it equal to zero. As we can see, when the modes are not 
coupled, that is, when A = 0, the evolution matrix reduces to the simple diagonal expression 

e iu 

[tA=o] = M = I I (22) 



e 



-in 



The phase difference between the uncoupled modes E + and E is equal to 2u. When A ^ 0, 
the evolution of the modes are coupled and the [U] matrix can be factorized as 



[U] = [U x =o}[M] = [u\[M\ (23a) 

\ / cos A - sin A e~ i{u - v) 
. A e i{ - u ~^ cos A 



(23b) 



Such a factorization will help us to estimate the TZ + and components of the different 
matrices. By using Eqs. fllOaj) . (11 Pel) and (TTTj) . we find 

K = 0, T+ = e-™, Q u = e~ 2iu , (24) 



so that 



and 



Hence, 




u ] — TFZ I , (25) 



K+ =tanAe i ("-"), T+ = -, M = 1. (26) 

COS /\ 



t m \n* M i 



Factorizing the free evolution phases as we did in Eq. ff23l) will allow us to point out a new 
phase expressed with the help of the © composition law. For this purpose consider the [U21] 
matrix which is the product of two unitary matrices \U\] and [U 2 ], 

[U 2 i] = [U 2 ][U 1 ]. (28) 

The factorization of the free evolution phases gives 

[U21] = [« 2 ][M 2 ][ui][Mi] (29a) 

= NfuilN^NtuiKM!] (29b) 

= [u 2 + u 1 \[M 2 {u 1 )\[M 1 \ (29c) 

= [m 2 + mi][M 2 i], (29d) 

where we have defined the diagonal matrix [u 2 + Ui] = [u 2 ] [u\] and noted that [M 21 ] = 
[M 2 (mi)][Mi] with 

( cosA 2 -sinA 2 e~ i(2lll+M2 ^ 2 A 
[M 2 {u l )] = [u- 1 %M 2 ][u 1 ]= , 9 * • (30) 

\sm\ 2 e< 2ui+U2 - V2) cosA 2 / 

If we use the definition ( fTTl) . we easily find 

©Mi = @Af 2 = ©M 2 («i) = 1- (31) 

^From the composition law Eq. (TX5T) we obtain 

©Mi©M 2 («i) +'^M 2 (« 1 )^M 1 1 + ^Mi'^M 2 (w 1 ) 

"Mai — =+ ~ — + — + , {M) 

1 + T^m 2 {ui) ^-Mi 1 + ^"Mi ^M 2 («i) 

or using the composition law definition in Eq. ([I]) 

©m 21 = K f^® n+M \ (33) 

i»21 "7?+ m"T? + v 1 

7t Mi w /c m 2 («i) 

It is interesting to note that m 2 i comes from the non-commutativity of the composition 
law ©. Although distinct, the two composite quantities TZ^h © ^m 2 («i) anc ^ ^m 2 (ui) ® "^Mi 
have the same modulus, so that 6m 21 is a pure phase term 

6 M21 = e" 2 ^. (34) 
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Finally, the whole phase term Qu 2 i = e 2l ^ 21 associated with the [C/21] matrix is 

9^ = e- 2 ^ = Q U1+U2 Q M21 = QM 21 e- 2t(ui+u *\ (35) 

which gives the phase 

02i = ui + u 2 + 4>. (36) 

The non-commutativity of the © law implies G m 21 7^ 1 in Eq. (1331) and is responsible for the 
additional phase appearing in Eq. (1361) . 

B. Examples of the physical meaning of the variable 

In the following we illustrate the meaning of the phase term by three examples from 
different fields of physics. 

1. Special relativity 

We first choose the composition of two non-parallel velocities v\ and v%- In this case the 
four elements A, B, C, and D of the matrix are respectively cosh(aj/2), sinh(aj/2) e~ tai , 
sinh(aj/2) e loli , and coshaj/2 where dj and Vi = tanha^ are respectively the rapidity and the 
velocity of the reference frame % for a given observer. Equations fllOap and (11 3 aft then give 
Vi = Kj = tanh (en/ 2) e iai and Vi = = tanh (oj/2) e~ m \ Here the phase a, gives the 
orientation of Ui with respect to an arbitrary axis belonging to the plane defined by the 
vectors vt and V2 in the observer reference frame. Because 0j = 1, Eqs. (1TB]) . ( 1321) . and ( 1331) 
give 

1 j_ 1? + -7?+ 1 _u T7 t/ 1 + tanh — e~ iai tanh — e ic * 2 
Q = 1 + K x Tt 2 = 1 + V1V2 = 2 2 ( 37 n 

21 1 + TZ+nt l + VV 2 1 + tanh ^ e iQl tanh ^ e— 2 ' 
1 2 2 2 

This expression is a pure phase term and can be written as 

e 21 = e~ 2i \ (38) 

where 2(f) is the Wigner's angle associated with the Thomas precession. Note that from 
Eq. (1371) we directly obtain the value of the Wigner angle. The real part of 621 gives 
immediately the known resultr^ 

/ ^1 i ^2 \2 / ^1 ^2 \9 

(1 + tanh — tanh — cos a) — (tanh — tanh — sin a) 

cos2(f) = ~, rk — rl , 2 ; : I — ~k ( 39 ) 

(1 + tanh — tanh — cos a) + (tanh — tanh — sin a) 
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where a = a% — a\. 

As was shown at the end of Sec. IIV A\ the Wigner angle comes from the non- 
commutativity of the composition law which mimics the non-commutativity of the Lorentz 
boosts. The © law can be easily used for the composition of any number of co-planar 
velocities. For example, for three referential frames, we obtain by iterating Eq. (TjQ) 

y + y + y + yy y 
W = V 3 © (V 2 © Vi) = _ 2 _ 1 _ . (40 

2. The optics of stratified media 

This example is from the optics of stratified media. If r$ and ti are the Fresnel reflection 
and transmission coefficients of the interface i (i — 1,2), and fa is the phase shift associated 
with the propagation of light, the four elements A, B, C, and D of the matrix (Q are 
respectively 1/U, (ri/U) e~ l<t>i , (fi/ti) e 1 ^*, and l/t». Equation ( llOah gives TZf = rie 1 ^, so 
that Eq. (I16al) gives the overall reflection coefficient^ of the two interfaces ([6]). In this 
case a phase term^ also appears which is strictly similar to the Wigner angle in special 
relativity. Its origin comes also from the non-commutativity of the © law which is related 
to the non-invariance of the problem when the two interfaces are exchanged. 



3. Light wave polarization 

As our last example, we consider two non-perfect polarizers. As explained in Sec. [TTJ the 
quality of the polarizer resulting from using successively two polarizers Pi and P 2 is given^ 
by Eq. (CQ). It is interesting in this case to calculate the value of Q 2 \. As was explained, 
02i expresses the non-commutativity of the two quantities which are composed. In special 
relativity finding two different results when calculating the resulting velocity of v\ composed 
with v 2 and of v 2 composed with v\ might have been surprising. It is not the case with 
polarizers. It is well known that the final polarization of a light wave going through the 
polarizer Pi and then through the polarizer P 2 is not the same as the final polarization of 
the light wave going first through polarizer P 2 and then through polarizer P\. We consider 
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explicitly the non-commutativity of polarizers. From Eqs. (J7|) and (fl8l) we obtain 

1 + tanh ^ e iai tanh ^ e" ia2 

e- 2in = V— ( 41 ) 

1 + tanh — e" mi tanh — e 4 " 2 
2 2 

In Eq. ( 14T1) 2Q is the angle between the polarization of light E 12 when going through the two 
polarizers in the order P\ and then P 2 and that of light E 21 when going through the polarizers 
in the order P 2 and Pi. For two perfect polarizers we expect to find 2Q = a 2 — ol\ = a- To 
verify this result, replace tanh — and tanh — by unity for perfect polarizers, and then the 
corresponding Eq. (l4Tj) for polarizers gives the expected result 

cos2fi = cos a. (42) 



V. DISCUSSION 



^From Eq. (1271) . we observe that there are redundancies of information in 2 x 2 unitary 
matrices. All information is contained in the first (or the second) column. Because of this 
redundancy, it is easy to understand why the use of the composition law (CQ) is easier and more 
rapid than using matrix methods such as transfer matrices. Moreover, as shown in Ref. 4], 
calculations converge more rapidly when the composition law is used. This rapid convergence 
comes from the fact that the denominator of the composition law is a normalization factor. 
Another useful aspect of the composition law is that it can be easily iterated 

K,...,i =K® (K-i © • • • W © Tit))- (43) 

As mentioned, this property leads to efficient algorithms for many kinds of problems. Also, 
Eq. ( f43|) is so simple that its analytic value can be directly given without any matrix calcu- 
lations. Equation ( 1431) is a complex generalization of the elementary symmetric functions of 
the mathematical theory of polynomials:— the numerator of Tl^... 1 is constituted by all the 
possible odd ordered products of the different 7Z^ factors, such that in each product, the 
1Z + and TZ + factors appear alternatively, the first factor always being 1Z + . The denominator 
of Tl^.„ lis constituted by all the possible even ordered products of TZf, such that in each 
product, the TZ + and TZ + factors appear alternatively, the first always being TZ + . If we limit 
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ourselves to two iterations, the value of TZ% ... x is directly given by 



(44a) 
(44b) 



1 + TZ[ K+ + K[ Kl, + 1l 2 n t 
Such a result can be useful for the case of S-matrices because the generalization of Eq. (J44 
allows us to write the S-matrix simply as 



S 



(Jit © Tit 



T- 



(T 2 - ® T{ 



A 



(T+®T+) Tl-®---{Tl 2 ®n^)j 



(45) 



It is well known that a number of physical processes are more adequately described by S'- 
matrices than by T-matrices (transfer matrices). Unfortunately, whereas T-matrices must be 
successively multiplied together, [M^...^] = [M n ] [M n _i] ■ ■ • [M 2 ] [Mi], such is not the case 
with S'-matrices. The composition law is consequently useful for S-matrices because our 
results show how to directly calculate the four elements of the overall S-matrix by iterating 
Eq. (05). 



VI. CONCLUSION 



The composition law of velocities in special relativity appears to be the natural way to 
add velocities that are subject to the condition \v\ < c. Its generalization in the complex 
plane leads to simple calculations of bounded quantities which would be otherwise difficult 
to calculate. We have shown how, for example, the Wigner angle in special relativity, the 
overall reflection coefficient of any multilayer, and the effect of any number of polarizers 
can be directly obtained from this general composition law. Also, we have shown that the 
generalization of the Einstein's composition law provides a natural way to compose scattering 
matrices. 
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